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Abstract 

The classic density-functional theory (DFT) formalism introduced by Hohenberg, Kohn, and 
Sham in the mid-1960s, is based upon the idea that the complicated A^-electron wavefunction 
can be replaced with the mathematically simpler 1-electron charge density in electronic struc- 
ture calculations of the ground stationary state. As such, ordinary DFT is neither able to treat 
time-dependent (TD) problems nor describe excited electronic states. In 1984, Runge and Gross 
proved a theorem making TD-DFT formally exact. Information about electronic excited states 
may be obtained from this theory through the linear response (LR) theory formalism. Begin- 
ning in the mid-1990s, LR-TD-DFT became increasingly popular for calculating absorption and 
other spectra of medium- and large-sized molecules. Its ease of use and relatively good accuracy 
has now brought LR-TD-DFT to the forefront for this type of application. As the number and 
the diversity of applications of TD-DFT has grown, so too has grown our understanding of 
the strengths and weaknesses of the approximate functionals commonly used for TD-DFT. The 
objective of this article is to continue where a previous review of TD-DFT in this series [Annu. 
Rev. Phys. Chem. 55: 427 (2004)] left off and highlight some of the problems and solutions from 
the point of view of applied physical chemistry. Since doubly-excited states have a particularly 
important role to play in bond dissociation and formation in both thermal and photochemistry, 
particular emphasis will be placed upon the problem of going beyond or around the TD-DFT 
adiabatic approximation which limits TD-DFT calculations to nominally singly-excited states. 

Posted with permission from the Annual Review of Physical Chemistry, Volume 63 ©2012 by 
Annual Reviews, http://www.annualreviews.org. 



CONTENTS 
INTRODUCTION 



WAVE-FUNCTION THEORY i 

Hartree-Fock Approximation 

Excited States 

Stability Analysis and Bond Breaking 

CLASSIC DENSITY-FUNCTIONAL THEORY 



Annu. Rev. Phys. Chem. 2011 1 1056-8700/97/0610-00 



Formalism . . . 
Approximations 
Scale-up .... 



TIME-DEPENDENT DENSITY-FUNCTIONAL THEORY 

Classic TD-DFT 

Applications 

Problems, Problem Detection, and Solutions 



GOING BEYOND THE ADIABATIC APPROXIMATION 

Current and Lagrangian Density-Functional Theory 

Optimized Effective Potential Approaches 

Time- Dependent Reduced Density Matrix Functional Theory . 

GOING AROUND THE ADIABATIC APPROXIMATION 

Reconciling TD-DFT and ASCF 

Spm-Flip TD-DFT 



29 

PERSPECTIVES [SC 



1 INTRODUCTION 



Space and spin coordinates 
have been denoted by, 
1,2,-- - = riCJi,r20-2,---- 

DFT: density- functional 
theory. A functional is a 
function of a function. It 
is designated by a square 
bracket notation. Thus 
u[/?](r) is simultaneously a 
functional of the function 
p and a function of the 
position, r. 

Unless explicitly present, 
we take fi = me = e = 1 
(i.e., we will use Hartree 
atomic units) throughout 
this article. 

TD: time dependent 
UV-Vis: ultraviolet-visible 
NLO: nonlinear optics 



ab initio: Latin meaning 
"from the beginning" used 
in quantum chemistry 
for first principles wave 
function calculations but 
often includes DFT in 
other contexts 
LR: linear response 



Even after the Born-Oppenheimer separation, the time-independent Schrodinger 
equation for N electrons in an "external potential" consisting of M nuclei and 
any applied electric fields, H^{1,2, ■ ■ ■ , N) = E^{1,2, ■ ■ ■ , N), is stih notori- 
ously difficult to solve. Yet such amazing progress has been made in solving this 
equation for increasingly complex systems since its introduction early in that 
century that we may well refer to the 20th century as the "Century of Quantum 
Mechanics." Of course, computers have helped quite a bit, but improved approx- 
imation methods based upon progress in our understanding of the mathematical 
and physical properties of the underlying objects of quantum mechanics has also 
been instrumental in finding more efficient solutions. One of the landmarks has 
been Hohenberg-Kohn-Sham ([UEj) density-functional theory (DFT), in which it 
is recognized that the complicated many-electron wave function, ^(1,2, • • • ,N), 
may be replaced with functionals of a simpler object — namely the charge density, 
p(ri) = iVEai/ /•••/|^(rio-i,2,--- ,iV)pd2d3---(iiV— in calculating proper- 
ties of the ground-stationary state. DFT was a mighty advance but it was to 
remain pretty much in the shadow of wave-function theory as far as quantum 
chemistry was concerned until into the 1980s when the tide turned. Approx- 
imate exchange-correlation functionals had advanced to the point where they 
could often compete with accurate ab initio wave-function methods, particularly 
for larger molecules of experimental interest. Nevertheless the restriction to the 
ground-stationary state was clearly disappointing for such popular applications as 
UV-Vis spectroscopy, nonlinear optics (NLO), and photochemistry. But, as with 
everything else, it seemed that "it was all just a matter of time." The formal basis 
of modern time-dependent (TD) DFT is the Runge-Gross theorem showing that, 
for a system initially in its ground stationary state exposed to a time-dependent 
perturbation, the time-dependent charge density, p{r,t), determines the time- 
dependent external potential up to an additive function of time ([3]). Once again 
there is a formally-exact functional which must be approximated but which then 
gives immediate access to NLO properties and, via linear-response (LR) theory, 
to excited-state information. The LR-TD-DFT equations introduced by one of 
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us (jU have been programmed in most quantum chemistry and quantum physics 
electronic structure codes and are probably currently the mostly widely-used 
method to treat the excited states of medium- to large-sized molecules. Many 
review articles have been written on various aspects of TD-DFT, some special 
volumes devoted to TD-DFT have appeared or will soon appear, and three books 
have been published or will soon be published (see the references cited in Ref. (jl]) 
as well as the newer Refs. dSJlT]).) The present review may be regarded as a follow 
up to an earlier review of TD-DFT in the present series ([8}). 

We take the point of view of theoretical physical chemists interested in spec- 
troscopy and photoprocesses and particular emphasis will be given to advances 
made over the last seven years. Figure [1] provides a schematic summary, sam- 
pling the breadth of chemical and physical phenomena that physical chemists 
might wish to describe. We will often want to refer to various features on this 
diagram. We will also often take pains to first review wave-function theory be- 
fore discussing (TD-)DFT. Our idea in doing this is to concentrate more fully 
afterwards on what is different in (TD-)DFT when compared to wave-function 
theory. 

This review is organized as follows. The next section reviews elementary wave- 
function theory. Section [3] reviews conventional time-independent ground state 
DFT. This is necessary in order to keep the chapter somewhat self-contained 
because TD-DFT shares many features in common with DFT and problems in 
approximate functionals used in TD-DFT may often be traced back to problems 
inherited from approximate functionals in TD-DFT. Section H] then introduces 
formal TD-DFT, the adiabatic approximation (A A) which defines conventional 
TD-DFT, typical applications, and then discusses problems with conventional 
TD-DFT, how to test for problems, and how problems can be fixed or at least 
attenuated. Some of these problems come from the very nature of the TD-DFT 
AA. TD-DFT should be exact provided the exact frequency- and initial state- 
dependent exchange-correlation (xc) functional of TD-DFT is used. This is not 
the same as the xc-functional of ordinary DFT which is simply borrowed by TD- 
DFT whenever the AA is used and so a true TD-DFT xc-functional should go 
beyond the AA to include frequency and possibly also initial state dependence. 
Section [5] discusses what is known about the exact xc-functional of TD-DFT 
and discusses attempts at creating and using approximate functionals which go 
beyond the TD-DFT AA. As these attempts are still in a more or less embryonic 
stage, it is also interesting to consider practical working methods to overcome the 
limitations of the TD-DFT AA without introducing any frequency dependence. 
This is the topic addressed in Sec. [6l Section [7] sums up with some perspectives 
regarding future developments in TD-DFT. 



Ref. (j8]): Marques MAL, 
Gross EKU, Annu. Rev. 
Phys. Chem. 55:427 (2004), 
"Time-Dependent Density- 
Functional Theory". 



AA: adiabatic approxima- 
tion 



xc: exchange-correlation 



2 WAVE-FUNCTION THEORY 



Modern DFT and TD-DFT resemble wave function theory in its simplest incar- 
nations. For this reason it is useful to review basic wave- function theory if only 
to establish some notation and to make sure that key concepts are fresh in the 
mind of the reader. 
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2.1 Hartree-Fock Approximation 

According to basic quantum mechanics, it is the wavefunction, ^, which de- 
scribes the fundamental state of the system and so we should aim to solve the 
Schrodinger equation. Except for exceptionally simple systems, we can only solve 
the Schrodinger equation approximately with the help of powerful approximation 
techniques such as the well-known variational principle. 

HF: Hartree-Fock. The Hartree-Fock (HF) approximation consists of using a trial wavefunction 
in the form of a single Slater determinant of orthonormal orbitals. Varational 

minimization of the energy subject to the orbital orthonormality constraint leads 

to the HF equation, 

/>,(!) = £^^^,(1), (1) 

given in a spin-orbital notation. Here the one-electron Fock operator, /, is the 
sum of the one-electron kinetic energy operator, t, the external potential, Uext(r), 
which is a multiplicative operator describing the interaction of the electron with 

the electric fields generated by the M nuclei and any applied electronic fields. The 

SCF: Self-consistent field. Fock operator also includes a self-consistent field (SCF) composed of a Hartree 
term, fjf(ri) = / /o(r2)/ri2 dr2, (also known as the Coulomb term), and an ex- 
change operator, Sj; (designated here as the exchange self-energy using the ter- 
minology and notation of Green's function theory.) As described in elementary 
quantum chemistry textbooks, 

±.m = - I ^^<A(2)(i2, (2) 

is a complicated integral operator, where 7(1,2) = X^j V'«(l)'^iV'i*(2) is the one- 
electron reduced density matrix (1-RDM) and the rij are orbital occupation num- 

bers. For completeness, we note the HF energy expression for the total energy, 

1-RDM: one-electron re- » 

duced density matrix. ^ = X! '^i{4>i\i\4>i) + / Vext{'r)p{r) dr + Eh[p] + E^Ij] , (3) 



where i is the one-electron kinetic energy operator, Vext is any potential external 
(i.e., the nuclear attraction and any applied electric fields) external to the system 
of A^-electrons, the Hartree energy, 

E»W = ^//«2^<;r,*„ (4) 

and the exchange energy, 

EM = -yP^^U2. (5) 

According to Koopmans' theorem ([9]), the occupied orbital energies, e^, in Eq. ([1]) 
may be interpretted as minus ionization potentials and the unoccupied orbital 
energies as minus electron affinities. That is, occupied orbitals "see" A'^ — 1 
electrons while unoccupied orbitals "see" A'' electrons. 

2.2 Excited States 

Let us now turn to the excited-state problem. We will begin with the equation-of- 
motion (EOM) formalism. This has the advantage of conceptual simplicity and. 
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though we will not show it here, can be directly linked to linear response theory 
p^ . The operator, |/)(0|, which destroys the ground state (|0)) and creates the 
Jth excited state (|/)) and its adjoint, |0)(/|, are solutions of the EOM, 



Evidently u = Ej—Eq is an excitation energy in the case of |/)(0| and oj 
is a de-excitation in the case of |0)(I|. 

We will now seek a solution of the EOM of the form, 



En 



(6) 

-Ej 

(7) 



EOM: equation-of-motion. 
npmh: n-particle/m-hole. 



That is we develop in a basis set of one-particle/one-hole (Iplh) excitation 
and de-excitation operators. This operator basis and the EOM operator "metric" 



defined by {A,B) = {^hf\ 
of the EOM, namely, 



\^hf) then allows us to develop a matrix form 



(8) 



A 


B 






1 




B* 


A* 


if] 


I- 


0-1 





where the notation X and Y is deliberately chosen to indicate that these objects 
are to be treated as column vectors in Eq. ([8]). The matrices are defined by. 



Aiajb 

where the coupling matrix. 



^i,j^a,b i^a 



pq,rs 



{PQlfnlsr) - {pr\fH\sq) 



(9) 



(10) 



Note that here and hence forth we adapt the molecular orbital (MO) index con- 
vention, 

a,b,c,- ■ ■ ,g,hi, j, k,l,m,n o,p,q, ■ ■ ■ ,y, z , (11) 

^^^^^^^v^^^^^^^ ^^^^^"v^^^^^^ ^^^^^^"v^^^^^^^ 

unoccupied occupied free 

where "free" means "free to be either occupied or unoccupied." We are also 
using a "lazy version" of second-quantized notation where p = hp and = a^, 
and Mulliken "charge cloud" notation. 



{pq\f\sr) 



V^;(l)Vg(l)/(l,2)V:(2)V'r(2) dld2, 



(12) 



where we have taken / to be the Hartree kernel, fni^, 2) = l/ri2. 

Equation ([8|) is often called the random phase approximation (RPA) for purely 
historical reasons. Some caution is in order when using this term since the 
RPA was originally developed in nuclear and solid-state physics without the sec- 
ond integral in the coupling matrix [Eq. (jlOp ]. Although we have chosen an 
EOM derivation of Eq. (|10p . this equation also emerges from LR-TD-HF and so 
may be referred to by that name (see for example, Ref. (jlOp pp. 122-134, 144- 
151). It is a pseudoeigenvalue problem with an operator "overlap matrix" on 
the right-hand side which is simplectic rather than positive definite as would be 
a true overlap matrix. It also has paired excitation, {X,Y), and de-excitation 



MO: molecular orbital. 



RPA: random phase ap- 
proximation. Sometimes 
one talks of RPA with ex- 
change (RPAE or RPAX) 
in order to emphasize 
inclusion of the second 
integral in Eq. (jlOp . How- 
ever it is often left up to 
the reader to figure out 
whether this integral is 
included or not in the term 
RPA. 
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TRK: Thomas- Reiche- 
Kuhn /-sum rule. 



SI: Systeme International. 



TDA: Tamm-DancofF ap- 
proximation. 

CIS: configuration interac- 
tion singles. 



TOTEM: two-electron 
two-orbital model 



MSM: multiplet sum 
method 



solutions, {Y*,X*), whose corresponding excitation/de-excitation energies, u, 
differ only by a sign change. The normalization of the solutions is given by 
= X'^X — Y^Y, which is typically positive for excitations and negative 
for de-excitations. Potential energy surfaces (PESs) may be calculated from the 
excited-state energies, Ej = Ehf + f^/- Since analytic derivatives are available 
for LR-TD-HF (jlip . automatic searching for critical points on excited-state PESs 
(Fig. [1]) is possible in this model. 

An important point that will only be mentioned here is that the LR-TD-HF 
equation provides not only excitation energies, loj, but also the corresponding 
oscillator strengths, 



2mf 
"3^ 



<^imr\I)\' 



(13) 



which may be calculated from the coefficient vector, {X,Y), and the dipole ma- 
trix. Oscillator strengths are pure numbers and, in a complete basis set (or other 
basis set satisfying the relation = +ih, satisfy the Thomas- Reiche-Kuhn 

(TRK) /-sum rule, X)/ // = -^j where N is the number of electrons. Equa- 
tion has been deliberately expressed in Gaussian, rather than atomic, units 
so that it may be more easily related to the molar extinction coefficient, e, in an 
absorption experiment. In SI units, the frequency spectrum is given to a first 
approximation by. 



6(^) 



^ 4m,cln(10)6o'^^^ " 

(6.94 X 10+is L/cm.s) ^5 

I 



(14) 



where Na is Avogadro's number and S{i') is a spectral shape function (typi- 
cally a Gaussian whose area is normalized to unity and whose full-width-at-half- 
maximum is determined from experiment.) 

An approximation which is frequently made is the Tamm-Dancoff approxima- 
tion (TDA), which consists of neglecting the B matrix in the LR-TD-HF equation, 
thus decoupling the excitations from the de-excitations. The TRK /-sum rule is 
lost when the TDA is made, but the resultant equation is variational in the wave 
function sense. In particular, the minimal acceptable model for the excited states 
of a closed-shell molecule is to take the wave function as a linear combination 
of Iplh (singly-excited) HE configurations, ^f^^ = '}2i,a'^'^^^HFCia,i = 0\(^hf- 
Use of the variational principle to find the coefficients which minimize the energy 
then yields the equation, 

AXi = uJiXi, (15) 

where Cj has been renamed Xj and the notations A = H — Ehf and loj = 
Ej — Ehf have been introduced in order to make the correspondence with the 
TDA of LR-TD-HF clear. Because of the variational nature of the TDA, for 
some applications, this "approximation" sometimes gives better results than the 
original RPA equation {vide infra). 

It is useful for interpretational purposes to restrict ourselves to the case where 
an electron is excited from an initial MO i to a target MO a neglecting any 
orbital relaxation. This two-orbital two-electron model (TOTEM) corresponds 
to the spin-preserving excitations in Fig. [2] and yields one singlet and one triplet 
excitation energy, as shown in Table [TJ 
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The CIS expressions of Table [T] may also be derived in a more direct way 
using the multiplet sum method (MSM) as fohows: The wave functions in Fig. [2] 
are ah eigenfunctions of the spin operator Sz but not necessarily of S'^. Linear 
combinations, ^s,Aisy needed to form simultaneous eigenvectors of the two 
operators, 

v&i,o = -^{<Pts-<P':i) , ^i,+i=<^tp, (16) 

where ^ifl, and ^i^+i are degenerate triplet states and $o,o is an open- 

shell singlet state. The triplet energy is easily expressed as the expectation value 
of a single determinant, 

Et = = |^|1>f ) = E[<^>^^] . (17) 

The MSM allows the singlet energy, Es, also to be expressed in terms of expec- 
tation values of single determinant wave functions. To do so, we note that, 

Es = (^o,oi^i^o,o) = ('5ri^i^r) + (^?ai^i<> 

Et = (^i,o|i?|^i,o) = ('5?a"|i^l$?.")-(«'?.1^l<). (18) 
Eliminating the cross term and using Eq. (jl7p gives us finally, 
Es = 2($-|i^|$-) - (cDf l^l^f ) 

= 2EMS]-E[^t], (19) 

where the last line emphasizes that the singlet energy has now been expressed 
uniquely in terms of the energies of single Slater determinants. Subtracting 
the ground-state HF energy, Ehf, and explicitly evaluating the various single- 

determinantal energies then leads to the CIS formulae in Table [TJ 

Unfortunately the CIS formulae in Table [1] turn out to be rather bad approx- STEX: static exchange 

imations for the singlet and triplet excitation energies because, as Koopmans' 

theorem tells us, HE orbital energies are preprepared to describe ionization and 
electron attachement rather than neutral excitations. It would be better to begin 
with a zero-order picture where both the initial and target orbitals "see" — 1 
electrons. This is the idea behind the static exchange method (STEX) (jlSp where 
the excitations of a neutral molecule that begin with MO i are carried out by first 
performing a calculation on the cation with one electron removed from the MO 
i (denoted by the abbreviated configuration, The associated unoccupied 

cation MO creation operators, a^{i~^), are then used together with the occupied 
neutral MO destruction operators, i, to carry out a variational calculation with 
a wave function of the form = J2aO'H^~^)'''^HF- Applying this idea to the 
TOTEM gives the STEX formulae, 

DM = Aaii~^) - li 

W5 = LOM + {ia\fH\ai) 

ujT = LOM - {ia\fH\ai) , (20) 

where the singlet, us, and triplet, wt, excitation energies have been expressed in 
terms of an excitation to a fictitious single-determinantal state of mixed symmetry 
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whose excitation energy, um, is equal to the difference between minus the electron 
affinity of the cation, Aa{i~^) and minus the ionization potential of the MO 
i, !{. Specific CIS values are given in Table [Ij 

It is interesting to ask what happens for a charge-transfer excitation between 
two neutral molecules located at a distance R apart, where R is assumed to 
be large. According to Eq. (j20p . we have ojt ^ ^ when the differential 
overlap ^j(r)^a(r) becomes negligeable. It is easy to see from Koopmans' theorem 
and the CIS results in Tabled that ujm ~ Aa — — R as would be expected on 
physical grounds. That is, the excitation energy required to tranfer an electron 
from MO i to MO a over a distance R is just the energy needed to ionize the 
electron from MO i plus the energy recuperated from adding it into MO a and 
from the electrostatic attraction of the two ions. 



2.3 Stability Analysis and Bond Breaking 



SODS: same-orbitals-for- 

different-spins 

DODS: different-orbitals- 

for-different-spins 



Most HF calculations are really restricted HF calculations with the same-orbitals- 
for-different-spins (SODS, as opposed to DODS, which stands for different-orbitals- 
for-different-spins.) In addition, they assume real MOs, each of which belongs 
to an irreducible representation of the molecular point group. Unrestricted HF 
consists of dropping the restriction of same orbitals for different spin, but still 
typically uses real MOs. As unrestricted HF calculations are more general than 
restricted HF calculations, cases arise where unrestricted HF calculations give a 
lower energy than restricted HF calculations by allowing spatial symmetry break- 
ing. Since the iterations in the HF SCF calculation typically preserve symmetry 
for closed-shell molecules, finding unrestricted HF calculations with lower energy 
than restricted HF calculations typically requires explicit symmetry breaking to 
be introduced during the iterations. On the other hand, such symmetry breaking 
can, for example, be a simple way to dissociate H2 into neutral atoms, Hf + H^, 
which is often a good approximation to the correct dissociation into a structure, 
[Ht + H| ^ H| + Hf], with equal probability for finding each spin on each atom. 
Generalized HF removes all restrictions on the HF orbitals. It turns out that 
there is an intimate relationship between restricted HF stability analysis and 
LR-TD-HF. 

Symmetry breaking will occur when an orbital unitary transformation, "(/"^(r) = 
e^^* ipri''^) = e*'^'^^"''*^'''i/'r(r) leads to energy lowering. 



Ex = Eo + — i W* W 



A B 
B A 



U 

u* 



R^ (A-B)i? + /^(A + B)/1 +0(3)(A) 



(21) 



where the A and B matrices are the same as previously defined in Eq. ([9]) except 
that they are real because we are beginning with a restricted HF solution with real 
orbitals. As Eq. ([8]) can always be put into the form (A + B) (A — B) Z = ujZ, 
where Z = X — Y, it is immediately evident that an imaginary LR-TD-HF 
excitation energy means that one or both of the matrices A -|- B and A — B have 
negative eigenvalues and thus that the restricted HF solution is variationally 
unstable because there is some choice of the vectors R and / which will further 
lower the restricted HF energy. The most common of these is the so-called "triplet 
instability" where the imaginary energy is associated with a triplet excitation 
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energy. However other types of instabilities are possible. For a more complete 
analysis see for example Ref. ()16p . 

Note that imaginary excitation energies can never occur when the LR-TD-HF 
equation is solved in the TDA. In fact, the variational nature of CIS does much 
to prevent the unphysical plunges of excitation energies which are often seen in 
LR-TD-HF as, say, singlet near instabilities associated with triplet instabilities. 
Thus the TDA is a partial solution to the LR-TD-HF instability problem. 

Physically LR-TD-HF is giving nonsensical excitation energies because it is 
based upon the response of an unphysical ground state wave function. In fact, 
while the HF approximation is often a good starting point for closed-shell molecules 
at their ground-state equilibrium geometry — i.e., it may be reasonable around 
the minima on the ground-state curve in Fig. [1] — the same approximation is 
notorious for failing at transition states where bonds are being made or broken. 
In this case, the minimum acceptable description is a wavefunction which is a 
linear combination of the reactant, and product, <I>p, trial wave functions, 
$ = ^rCh + ^rCr. Typically this involves some sort of configuration mix- 
ing between the ground-state determinant and a two-particle/two- hole (doubly) 
excited state. Though there is no transition state on the ground state curve, 
the dissociation of the H(ls)-H(ls) u-bond is a textbook case. The minimum 
acceptable correctly dissociating description of the ground state wave function, 
satisfying all the spatial and spin symmetry rules, is a trial function of the form, 
^ = \crgag\Ci + \cruO'u\C2, where (Jg is the familiar bonding MO and au is the 
familiar antibonding MO. The coefficients Ci = l/\/2 = —C2 give the proper neu- 
tral, H2 —7- [Ht + H| -H- H| + Ht], dissociation limit of the ground X^Tig state. 
The combination Ci = l/\/2 = C2 gives the proper ionic dissociation limit, 
H2 ^ [H - + H+ o H+ -FH:"], of the "two-particle/two-hole" ("2p2h") l^Sg 
state. Though the l^S^ state is dominated in this simple model by the 2p2h 
configuration at the ground-state equilibrium geometry, the fact that neither Ci 
nor C2 are zero in the dissociation limit indicates that configuration mixing of 
the ground and 2p2h excited-state configurations is occuring. In fact the simplest 
wave-function model for breaking of any homoleptic single bond involves mixing 
of the ground-state configuration with a 2p2h configuration. 

3 CLASSIC DENSITY-FUNCTIONAL THEORY 

It is difficult (possibly impossible) to speak of TD-DFT without first saying a 
word about ordinary DFT. This is the classic DFT of the ground stationary 
state. It should also be the static limit of TD-DFT. Since our focus in this 
chapter is not on DFT, which is already fairly well known, but on TD-DFT, we 
shall try to keep this section short. Readers interested in more information about 
DFT may choose to consult any one of a number of books on the subject such 

as, for example, that of Dreizler and Gross (|17p . 

Before going further, it is important to make the distinction between DFT DFA: Density-functional 
(density- functional theory) and density- functional approximations (DFAs). DFT approximation. 

is a formalism which provides existence theorems showing that there is something 

well defined to be approximated. However from a practical point of view, DFT 
is pretty useless unless practical DFAs can be found which are simpler to use 
than the wave-function theory of corresponding accuracy. The search for prag- 
matic approximations, ideally guided by rigorous formal DFT but occasionally 
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borrowing from wave-function theory, pervades all of density-functional culture. 



3.1 Formalism 

Hohenberg and Kolin set down the basic formalism of modern DFT in two theo- 
rems ([T]). The first Hohenberg-Kohn theorem states that the ground state charge 
density of a system of interacting electrons determines the external potential up 
to an arbitrary additive constant. Since integrating over the charge density gives 
the number of electrons, this implies that the ground-state charge density de- 
termines the Hamiltonian up to an additive constant (i.e., the arbitrary zero of 
the energy) and hence determines pretty much everything about the ground- and 
excited-states of the system. 

In fact, the ground state density must even determine the dynamic linear re- 
sponse of the ground stationary state, 



hOO 



X(ri,r2;a;)= / 6+^^'^^) Z^^^' . d(ti - t^) , (22) 

J-oo 0Vappl{r2,t2) 

since this can be expressed in terms of the ground- and excited-state wave func- 
tions and excitation energies, 

^^ (0|p(ri)|/)(/|p(r2)|0) (0|/)(r2)|/)(/|p(ri)|0) ^ 
X[ri, r2; ui) = 2^ [ — ■ — , 23) 

where p is the density operator and rj = 0^ is an infinitesimal, needed to guarantee 
causality: 6p{ri,ti)/6vappi{r2,t2) = if ^2 > ^i- 

Counterbalancing this marvelous theorem are the representability problems. 
That is, given a candidate ground-state charge density which integrates to N 
electrons: (i) How can we be sure that it comes from an A^-electron wave function 
(pure-state representability) or at least from an A^-electron ensemble density 
matrix (ensemble N representability)? (ii) And how can we be sure that it is 
the ground-state density for some external potential {v representabilty)? It turns 
out, with few exceptions, that most reasonable candidate charge densities are 
pure-state representable, but determining the u-representability of a candidate 
charge density is far from obvious. 

There is also the problem of how to actually calculate something nontrivial 
such as the ground-state energy from a candidate charge density. The second 
Hohenberg-Kohn theorem solves this problem by using the variational principle 
to show that, Eq < F[p] + J Vext{^)p{^) dr^ where the functional F[p\ is uni- 
versal in the sense that it is independent of the external potential. Equality is 
achieved for the ground-state charge density. In this case, we can get around 
the u-representability problem using the Levy-Lieb constrained search formal- 
ism, F[p\ = inf.i,_j.p(^'|T -|- ygel^)) which says to search over all wave functions 
(or ensemble density matrices) which produce the trial density. In this sense, 
the Hohenberg-Kohn functional is not unknown, but it is too complicated to be 
practical and so approximations must be sought. 

Such approximations fall under the category of orbital- free theories, of which 
Thomas-Fermi theory is the classic historical example. They tend not to be as 
accurate as approximations based upon the Kohn-Sham (re) formulation of DFT. 
The main problem in making a DFA for an orbital-free theory is in finding an 
approximation for the kinetic energy. The Kohn-Sham formulation of DFT gets 
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around this problem introducing a fictitious system of noninteracting electrons 
with the same charge density as the real system of interacting electrons. The 
total energy expression is exactly the same as that in the HF approximation 
[Eq. ([3])] except that the HF exchange energy, -E^M) is now replaced with the 
exchange-correlation (xc) functional, Exc[p], and the total electronic energy is 
exact when the xc-functional is exact (in contrast to HF which is always an ap- 
proximation). Note that the xc-functional, Exc[p] = F[p]— EH[p] — ^ni^^p{^\T\^) 
where $ is a single-determinantal wave function, contains not only the exchange 
and correlation energies familiar from wave function theory but also the difference 
of the kinetic energies of the interacting and interacting systems. Also, before 
going any further, let us follow modern practice and generalize the original Kohn- 
Sham DFT to spin DFT so that Exc = Exc[pa, P/b]- Minimizing the total energy 
with respect to the orthonormal orbitals of the noninteracting system then leads 
to the Kohn-Sham equation which differs only from the corresponding HF orbital 
equation ([1]) by the replacement of the HF exchange self-energy [Eq. ([2])] by the 
xc-potential, 

SExc[Pa,P/3] 



Vxcipa, P^]{r) 



^Pa{r) 



(24) 



It is now well established that there is a particle number derivative discontinuity 
(PNDD) in the exact xc-potential, whereby adding a fraction of an electron into 
a higher-energy orbital leads to a drastic change in the asymptotic behavior of 
the xc-potential and a nearly rigid shift of the xc-potential in the "bulk region" 
where the density is large. This may be explicitly verified within the optimized 
effective potential (OEP) model which is most easily understood in its classic 
exchange-only version where it is often called exact exchange (EXX). The EXX 
model finds the potential whose orbitals minimize the HF energy expression. 
Alternatively the OEP potential is that for which the linear response of the 
density is zero when replacing the HF exchange self-energy with a local potential. 
Thus the OEP provides an excellent approximation to the true Kohn-Sham x- 
potential, complete with PNDD. This is easiest to see within the Krieger-Li- 
lafrate (KLI) approximation to the EXX potential for which we will now give a 
heuristic derivation: Let us assume that the HF and KS orbitals are about the 
same on average. If they really were identical, then we could immediately write 
that, (Tjx — f^x(r)^ V'j(r) = (^i^^ ~ ^f""^) ^«(r)- However this is only supposed to 
be true on average so we should multiply by ni^^{r) and sum over the index i. 
After rearrangement, this leads to. 



j:im^*ir)txMr) 
p{r) 



+ 



HF 



p{r) 



(25) 



The first term on the right-hand-side is Slater's historic definition of the local 
exchange potential. The second term is the PNDD correction. It is there because 
the asymptotic behavior of the charge density depends upon the value of HOMO 
energy. This second term introduces a rigid shift to compensate for the change in 
the asymptotic behavior every time a new HOMO is occupied. Such discontinuous 
behavior is very difficult to incorporate into practical DFAs other than by an 
OEP-like approach. 

With the Kohn-Sham formulation, DFT also gains a new represent ability prob- 
lem, namely the problem of noninteracting u-representability (NVR). This is the 



PNDD: particle number 
derivative discontinuity 
OEP: optimized effective 
potential 

EXX: exact exchange 
KLI: Krieger-Li-Iafrate 



NVR: noninteracting 
representability. 
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requirement that there be a noninteracting system integer occupation numbers 
whose ground state has the same density as that of the interacting system. Vi- 

olation of NVR would show up as a "hole below the Fermi level" meaning that 

HOMO: Highest occupied the lowest unoccupied molecular orbital (LUMO) is lower in energy than the 
molecular orbital. highest occupied molecular orbital (HOMO) in violation of the assumed Fermi 
LUMO: Lowest unoccupied statistics for the ground state of the noninteracting system. Now suppose ev- 
molecular orbital. ery system were NVR. Then we should always have an exact spin-restricted 
(same-orbitals-for-different-spin) solution for dissociating molecules with a sin- 
glet ground state belonging to the totally-symmetric representation of the molec- 
ular symmetry group. That is, symmetry breaking should never occur for the 
exact functional. This is because the exact ground state spin a and spin /3 den- 
sities are identical, implying that v^^{r) = v^^{r), and hence that the orbitals 
can always be chosen identical for the different spins (llSp . It now seems likely 
that not every system is NVR (see Ref. ()19p and references therein). This leads 
naturally to fractional occupation numbers as is easily seen in the TOTEM. Con- 
sider a spin-wave instability, where the Kohn-Sham determinant takes the form, 
= — A^^j + iXtlja, \/l — A^-i/Sj + iXi^al, where the absence (presence) of 
an overbar indicates spin a (/3). The total energy may then be expanded as 
Ex = Eq + 2A^(ea — Cj) + 0{X^). The spin-wave instability will lower the en- 
ergy whenever the LUMO energy, e^, is lower than the HOMO energy, €a- (See 
Ref. ()20p for a similar example within the HF model.) However it is then easy to 
demonstrate that the orbitals have fractional occupation number by looking at 
the corresponding 1-RDM operators, 

% = (1 - X^)\i){i\ + iXVl - A2 {\a){i\ - \i){a\) + A>)(a| 
7^ = (1 - - iXVl - A2 {\a){i\ - \i){a\) + A>)(a| 

7 = 7a + 7/3 = 2(l-A2)|i)(i|+2A2|a)(a|. (26) 

These observations are consistent with the traditional assumption that an en- 
semble formalism with fractional occupation number is needed in cases where 
NVR fails. An important theorem in this context is that all fractionally-occupied 
orbitals must be degenerate with lower energy orbitals being fully occupied and 
higher-energy orbitals being completely vacant (see e.g., Ref. ()17p . pp. 55-56). 
Surprisingly there is some indication that the noninteracting system may become 
degenerate in such a theory for a range of molecular configurations where the 
interacting system is nondegenerate ()2ip . One is then possibly faced with the 
problem of minimizing the energy with respect to unitary transformations within 
the space spanned by the degenerate occupied orbitals with different occupation 
numbers. 



3.2 Approximations 



Pure density functional: 
One depending only upon 
the density, excluding any 
orbital dependence. 



DFT by itself is hardly useful without DFAs. John Perdew has organized func- 
tionals into a Jacob's Ladder (Fig. [3]). The basic idea is that gradually including 
more information about the local density, the reduced gradient, the local kinetic 
energy, etc. allows users (angels going up the ladder) the added flexibility needed 
to create more accurate functionals. Only the first two rungs (or three if the local 
kinetic-energy is excluded from the third rung) constitute pure density functionals 
— that is, functionals which stick to the original spirit of Hohenberg-Kohn-Sham 
DFT and consider only the density as the working variable. For this reason it 
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might be better to refer to the higher rungs as generahzed Kohn-Sham. It must 
be emphasized that, while greater accuracy is not guaranteed, there is a general 
tendency towards more accurate functionals when going to higher rungs of the 
ladder. However it is important that users also be able to descend the ladder 
when necessary as calculations necessarily become more complex and expensive 
on the higher rungs of the ladder. 

One of us has generalized Jacob's ladder to Jacob's jungle gym (i22]l by includ- 
ing additional ladders corresponding to higher xc-derivatives. For pure density 
functionals these would be defined by Eq. and. 



S Exc[Pa,Pl3] 

^P<7i(ri)5/9a2(r2) 



5---(ri,r2,r3) = /'f^^^"\f . (27) 

where the xc-kernel, fxc, and xc-hyperkernel, gxc, are needed for static response 
properties. There are other analogous quantities at each of the higher rungs. This 
type of representation is useful to keep in mind because derivative quantities are 
always harder to approximate than the quantities whose derivative is being taken. 
It may thus sometimes be appropriate to combine different approximations from 
different ladders when calculating certain types of quantities. 



3.3 Scale-up 

Since Kohn-Sham DFT resembles the HF approximation which is the basis of 
much of the development of (non-DFT) ab initio quantum chemistry, it is rela- 
tively easy to compare the two approaches and even to incorporate a Kohn-Sham 
solver into a wave-function ab initio approach. This makes especially good sense 
if DFAs are viewed as a way of extending the accuracy of wave-function ab initio 
approaches to larger systems then would otherwise be possible. In this sense, the 
problem of how to scale-up calculations to ever larger systems is highly relevant 
to the success of DFT and DFAs. However space limitations prevent us from 
saying more than only a few words about scale-up strategies. 

It has long been noted that the lower rungs of Jacob's ladder (i.e., the LDA 
and the GGAs) more closely resemble the Hartree than the HF approach. As 
such, rather than have a formally 0^^\N) scaling, the introduction of density- 
fitting functions could lead to a formal 0^^\N) scaling which could be further 
improved upon by integral prescreening. Beyond a certain system size (and hence 
once a certain level of computational power is reached), special algorithms can 
be used such as those which scale asymptotically as 0{N). These algorithms 
are especially easy to apply in the case of DFT and makes linear-scaling DFT 
extremely attractive for large-scale "chemistry on the computer" (j24p . 

A second scale-up strategy is the subsystem approaches where formal DFT 
(|25ll26p is extended to describe subsystems of larger systems — such as molecules 
hydrogen bonded to other molecules, molecules in solution, or even fragments 
of molecules covalently bonded to the whole — and corresponding DFAs are de- 
veloped. Not only does such an approach allow the treatment of larger systems 
by allowing a Kohn-Sham description of a subsystem to be embedded within a 
larger surrounding system described, say, by using an orbital-free DFA, but the 
subsystem approach is also attractive from an interpretational point of view as it 
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corresponds nicely, say, with the idea of functional groups in organic chemistry 
and active sites in biochemistry. An intriguing development is the successful 
implementation of subsystem density-functional theory for the case that the sub- 
system is covalently bonded to the larger system (j27y28p . 

DFTB: Density- functional A final strategy which is worth mentioning is the development of the density- 
tight-binding, functional tight-binding (DFTB) model with parameters obtained from Kohn- 

Sham calculations (|29p. While loosing some of the rigor of a true Kohn-Sham 

calculation, DFTB is a next logical step in a multiscale strategy to larger spatial 
and temporal scales than can normally be obtained by true Kohn-Sham calcula- 
tions. 



4 TIME-DEPENDENT DENSITY-FUNCTIONAL THEORY 



ASCF: Method consisting 
of calculating the energy 
difference of two SCF 
calculations. 



The classic Hohenberg-Kohn-Sham DFT is a theory of the ground stationary 
state. However, as mentioned in the introduction, chemistry is not just a mat- 
ter of the ground stationary state, but also of UV-Vis spectroscopy, NLO, pho- 
tochemistry, and other applications involving either time-dependent electronic 
properties and/or electronic excited states. 

In 1999, Singh and Deb reviewed the prospects for treating excited states within 
a DFT framework (I30p . Up to 1995 (j3]), the most popular method was the ASCF 
approach and its MSM variant. Since the classic Hohenberg-Kohn-Sham DFT 
can handle electronic ground states, it is possible to carryout self-consistent field 
(SCF) calculations for the ground states of a neutral molecule, its anion, and its 
cation. By taking energy differences, one has a rigorous method for calculating the 
first ionization potential and electron affinity. By adding or removing electrons 
from other orbitals, one has at least an ad hoc way to estimate other ionization 
potentials and electron affinities. For some of these ionization potentials and 
electron affinities, the method may even be considered to be reasonably rigorous 
to the extent that we can accept that the DFAs developed for the true ground 
states also should apply to the ground state of each symmetry, particularly if a 
single-determinant is a reasonable first approximation. Within the same assump- 
tion, the lowest triplet excitation energy of a closed-shell molecule may also be 
calculated reasonably rigorously within the ASCF method. For open-shell singlet 
excited states, the Ziegler, Rank, and Baerends (jl2p and Daul (j3Tl) adapted the 
MSM method as an excited-state DFA. The result in the TOTEM is given in 
Table [1] after linearization. The difference is that minus the ionization potential, 
Ij, of orbital i and ininus the Glcctron affinity, .A^ (i ), are given by expresions 
which are very different than those given by Koopmans' theorem. In fact, unlike 
the case of Hartree-Fock, one would expect Kohn-Sham orbitals to be prepre- 
pared to describe excitations in the sense that for pure DFAs, the occupied and 
unoccupied orbitals "see" the same potential and hence the same number of elec- 
trons. This is particularly clear in the expression for the mixed-state excitation 
energy, um, which is just a simple orbital energy difference when the densities 
of the initial (i) and final (a) orbitals are the same. The additional term is a 
correction for relaxation of the orbital densities during the excitation process. 
There is thus a lot of good physics in the MSM DFA for excited states. However 
the major difficulty with this method is that the form of the excited-state wave 
function must be estimated, typically based upon symmetry arguments. It thus 
only a first-order estimate in many cases and is also more difficult to automate 
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for excited-state calculations than is TD-DFT. 
4.1 Classic TD-DFT 

The question is to what extent the time-evolution of the charge density can be 
treated without having to solve the full the time-dependent Schrodinger equation, 

Hm{t) = z^*(t) , (28) 

for the A^-electron wave function, ^{t) = ^{1,2, N,t). Before discussing 
this question, it is useful to recall that the time-dependent problem is an initial 
value problem: Given ^'o = ^'(^0)1 propagate the wave function forward in time 
to obtain ^{t). Thus the wave function at time 4 is a functional of the wave 
function at time to, ^'[^'o](t). 

Since the time-dependent A^-electron Schrodinger equation can rarely be solved 
exactly, approximate solutions are often sought. This is typically done using the 
Frenkel-Dirac stationary action principle (|32 p 33 t[3^ . which states that making 
the action, 

A[^;^o]{t,to) = £ (^m^o]it%-^mmt')) - Eimt')) dt', (29) 

stationary with respect to variations of ^[^o](i') on the interval t' G ito,t) and 
subject to the constraints (5^'[^'o](to) = ^^I'l^oKO = 0, implies that the time- 
dependent Schrodinger equation is satisfied over the time interval {to,t). Here, 

E[^om = m^omiHM^om) , (so) 

is the "instantaneous energy" which can in principle also contain information 
about system dynamics from previous times. This emphasizes the main point of 
this stationary action principle, mainly that an approximate equation of motion 
may be developed from any approximate instantaneous energy expression. Vig- 
nale pointed out a problem with the original Prenkel-Dirac principle, namely that 
causality (i.e., ^'[^o](i) is determined by the boundary condition ^q) means that 
we are not in general free to set (^^'[*I'o](i) = 0. When this is taken into account 
the new Frenkel-Dirac- Vignale stationary action principle p5|) is that, 

6A[^; ^o]{t, to) = im^o]{t, to)m^o]{t, to)) , (31) 

subject only to the obvious constraint that 5^'[^'o](to) = implies that the time- 
dependent Schrodinger equation is satisfied over the time interval {to,t). Note 
that the right hand side of Eq. (j3ip set to zero in the original Frenkel-Dirac 
formulation. Fortunately the two stationary action principles often give the same 
result ()35p . which is good since the Vignale correction is rarely invoked. Note also 
that Lagrange multipliers may easily be added to apply the constraints frequently 
needed when deriving equations within a given model. 

Now let us return to the question of the extent to which the time-dependent 
wave function may be replaced by the time-dependent charge density. It is nat- 
ural to replace the instantaneous energy in the Frenkel-Dirac action with the 
Kohn-Sham energy and the wave function with the Kohn-Sham determinant. 
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The Frenkel-Dirac stationarity condition then leads to the time-dependent Kohn- 
Sham equation, 

-^V^ + Vext{t) + VH[p{t)] + Vxc[Pa{t),Pi3{t)] 

and this equation was used as early as 1980 ()36p to calculate photoabsorption 
cross-sections of rare gas atoms. However the formal justification did not come 
until later. 

Formal TD-DFT is usually traced back to the classic paper of Runge and Gross 
([3]) which tried to firm up earlier work on the same subject. Thus Runge-Gross 
TD-DFT is two decades younger than its stationary ground-state older sibling, 
the Hohenberg-Kohn-Sham theory. The Runge-Gross paper actually contains 
four theorems. The first is the Runge-Gross theorem which states that the time- 
dependent charge density, p{r,t), together with the initial wave function, 
determines the external potential up to an additive function of time. This implies 
that the wave function is determined up to an arbitrary phase factor, ^{t) = 
g-«0(*)\]/[^^ The proof proceeds in two steps. In the first step, it is shown that 
any Taylor-expandable time-dependent external potential is determined by ^'o 
and by the current density, j(r, f). In the second step, it is to show that p(r, t) 
is determined by j(r,t). A weakness in the second step of the proof was noted 
by Xu and Dhara (j37p but can be excluded for physically-reasonable external 
potentials which may be viewed as being made up of point charges (ISSpSOh . 
Although the TD-DFT community is very self- critical, the proof of this theorem 
has thus far withstood the test of time. The second theorem in the Runge-Gross 
paper proposes a hydrodynamical formulation of TD-DFT. The third theorem 
proposes the Frenkel-Dirac action can be made the basis for a variational theorem 
somewhat analogous to the second Hohenberg-Kohn theorem. Unfortunately this 
theorem is wrong as stated because (as is always the case with the Frenkel-Dirac 
principle) it does not handle causality correctly. Many subsequent reformulations 
of the action have been proposed including one by Rajagopal involving the Berry 
phase (|40p , formulations involving the Keldysh action (j41U42p , and Liouville space 
pathways ()43p . However the use of the Frenkel-Dirac- Vignale action (jSSp seems 
to be by far the simplest solution. The fourth theorem in the Runge-Gross 
paper proposes a Kohn-Sham formulation of TD-DFT. This is actually a very 
significant statement because it clarifies the nature of the time-dependent xc- 
potential. Whether written as the derivative of an xc-action, 

Vxc[Pa,Pfi;'^o,'^o]{r,t) = , (33) 

OPa{r,t) 

or simply postulated to exist, the xc-potential in time-dependent Kohn-Sham 
theory has a formal dependence on the initial wave functions of the real inter- 
acting, ^Q, and noninteracting Kohn-Sham, $0, systems of A^-electrons. This 
may be eliminated using the first Hohenberg-Kohn theorem if the initial state 
of the system is the ground stationary state. Even so, we do not recover the 
xc-potential of Eq. ()32p . where time is merely considered to be a parameter and 
the derivative is taken at fixed time. This is the important AA of TD-DFT which 
assumes that the xc-potential reacts instantaneously and without memory to any 
temporal change in the charge density. The adiabatic approximation is a local 
approximation in time. The exact TD-DFT xc-potential should have memory 
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in the sense of depending upon the past history of the charge density and is 
thus a functional of a function of not just (r, a) but also of t. The adiabatic 
approximation may be taken as the definition of conventional TD-DFT and in 
practical terms predated the Runge-Gross paper. The gauntlet thrown down by 
Runge and Gross to challenge us is to find how to include memory effects in a 
time- dependent xc-functional. 

The Runge-Gross formulation of TD-DFT was limited only to applied electric 
fields. This is sufficient for many problems because magnetic field effects are 
often significantly smaller than electric field effects in experimental studies of 
molecular properties, but it would be better to have a gauge-invariant theory 
treating electric and magnetic fields on the same footing. Ghosh and Dhara have 
extended TD-DFT to handle time-dependent magnetic fields by introducing a 
current-density dependence and an exchange-correlation vector potential ()44p . 
Since that time, TD current density-functional theory (CDFT) has developed 
considerably (see, e.g., Ref. (|45p for one recent paper.) Alternatively it is possible 
to formulate and carryout calculations within the framework of 4-component 
relativistic TD-DFT ([i6ll . 

4.2 Applications 

In the quarter century since the Runge-Gross paper ([3]) and the introduction 
of TD-DFT into quantum chemistry over a decade ago (jH), there have been a 
very large number of interesting applications of this theory. We will only try to 
highlight a few here that we find to be especially interesting. 
4.2.1 Absorption Spectra The major application of TD-DFT is for the 
calculation of electronic absorption spectra using linear response theory. Ear- 
lier (Sec. 13. ip we argued that the first Hohenberg-Kohn theorem guarantees that 
the dynamic linear response of a system of A^-electrons is determined by the 
ground-state charge density, but we had no practical way to construct this re- 
sponse. TD-DFT provides just such a constructive procedure through the dy- 
namic linear response of the charge density. The calculation consists of prop- 
agating the Kohn-Sham orbitals in the presence of a small dynamic perturba- 
tion, Vappi{r,t) = £{t) • r, and calculating the induced dipole moment, fJ-ind{t), 
which is the difference between the time-dependent and permanent dipole mo- 
ments of the system. The dynamic polarizability, a is defined by the linear term, 
fJ-indit) = J a{t — t')£{t). Applying the Fourier transform convolution theorem al- 
lows us finally calculate the dynamic polarizability as, a{uj) = fi^uj) / £ (oj) . Since 
the dynamic polarizability has the sum-over-states (SOS) form. 



where uj and // are respectively vertical excitation energies and the correspond- 
ing oscillator strength, the Laurentzian broadened stick spectrum can finally be 
obtained as, 



Such a procedure has been implemented, for example, in the Octopus code (j47p . 
and has the advantage of being able to go to calculate the absorption spectrum 
of a very large molecule molecule over a wide range of energies, albeit with only 
moderate spectral resolution. 



CDFT: current density- 
functional theory 



SOS: Sum-over-states. 




(34) 




(35) 



vr 



18 



Casida <fe Huix-RotUant 



Alternatively one of us dH) showed how LR-TD-DFT could be formulated 
through an analytic treatment of the poles of the dynamic polarizability to look 
like the RPA equation [Eq. ([8])] thus allowing a rapid implementation of TD-DFT 
in ab initio quantum chemistry wave function codes. Thus this is the method now 
found in nearly every internationally important quantum chemistry and quantum 
physics code. The derivation of Ref. ([4]) did not make the adiabatic approxima- 
tion and allowed for (fixed) fractional occupation numbers. The matrices A and 
B are (^-dependent in the general theory and the number of solutions exceeds the 
dimensionality of the matrix problem. However, for brevity, we will only write the 
matrix elements of the A and B matrices here in the adiabatic approximation, 

Aiaa,jbT = 5^,T5i,j5a,b{e.aa " eja) + {ia\fH + fxc\bj) 

Biaa,bjr = {ia\fH + f^^jb) ■ (36) 

Once LR-TD-DFT equations are formulated for TD-DFT, it is natural to use 
the TDA and this approximation was introduced into the TD-DFT literature by 
Hirata and Head-Gordon ()48p . Just as the variational nature of CIS is helpful 
in circumventing, or at least, reducing problems associated with triplet and near 
singlet instabilities in the HF RPA, so too the TDA is highly useful when calcu- 
lating TD-DFT excited-state potential energy surfaces (ll3p i8 t[T9]l . However one 
should be aware that the cost of the TDA is a loss of oscillator strength sum rules 
(j49]) with possible implications for the accuracy of corresponding spectra (f50|l . 

We would now like to say a few words about why TD-DFT works from a 
physical point of view by looking at the TOTEM. At least on the lower rungs 
of Jacob's ladder (LDA and GGAs) the Kohn-Sham orbitals are preprepared to 
describe electronic excitations since the occupied and unoccupied orbitals "see" 
the same xc-potential and hence the same number of electrons. This means 
that, except for charge transfer excitations which have large amounts of density 
relaxation. The TOTEM often applies fairly well in TD-DFT. The corresponding 
TDA TD-DFT singlet and triplet excitation energies are given in Table [TJ An 
examination of the size of the various integrals then shows that, 

tJT < ea - ej < , (37) 

with all three values becoming identical for Rydberg states where the differen- 
tial overlap '(/'i(r)V'a(r) becomes negligible. At lower energies we should expect 
to recover the STEX formulae. As shown in Table [H this is possible, at least 
roughly, within a sufficiently accurate model such as the optimized effective po- 
tential (OFF) model p4p . Direct calculation shows that the Koopmans' theorem 
interpretation of the occupied orbital energies as the negatives of ionization po- 
tentials holds better than it does for HF. Thus li = e^. However the energies 
of the unoccupied orbitals do not represent the negatives of electron affinities 
because the occupied orbitals "see" — 1 rather than N electrons. The quan- 
tity Aa{i~^) is minus the electron affinity of the orbital ipa for the cation with 
an electron removed from orbital ipi. The electron thus "sees" — 1 electrons 
and so Aa{i~^) is to a first approximation equal to ta of the A-electron neu- 
tral. However one Hartree and one exchange interaction with the electron in 
orbital tpi have been removed, so we now have to update the approximation to be 
Aa{i~^) = ea + (a^l/i? + /xc"N^)- (Note that the Hartree interaction is estimated 
in TD-DFT by the integral {ai\f^f\ia).) 
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4.2.2 Photochemistry The popularity of TD-DFT for modeling absorp- 
tion spectra is a strong incitation to extend TD-DFT to related photophenomena 
shown in Fig. [TJ We are immediately confronted with the fact that TD-DFT 
provides only excitation energies, coj = Ej + Eq, so that excited-state PESs 
must be calculated by adding the excitation energy obtained from TD-DFT 

T D — D FT 

to the ground state energy obtained by static ground-state DFT: Ej = 
Eq +ujj ~ . Exploring excited-state PESs is facilitated by the implemen- 
tation of analytic gradients for TD-DFT excited states, VEJ^~^^'^ = V-E^^^ + 
VlVj , in many electronic structure codes. 

The first and easiest of these phenomena that can be treated is fluorescence. 
The TD-DFT fluorescence energy is a TD-DFT singlet excitation energy calcu- 
lated at the optimized geometry of the singlet excited state. The corresponding 
Stokes shift is the difference between the absorption and fluorescence energies, 

(^STOKES = ^ABS — ^FLUO- 

More challenging is to use TD-DFT to investigate possible photochemical re- 
action mechanisms. While TD-DFT is now well-established as a part of the 
photochemical modelers' toolbox, applications to studying the entire course of 
a photochemical reaction are relatively new. The reader is referred to a recent 
review on the subject by Casida, Natarajan, and Deutsch (jSip . Suffice it to say 
here that it is now increasingly possible to go beyond Ehrenfest dynamics or 
the simple pathway method and carry out mixed TD-DFT/classical trajectory 
surface hopping calculations to get an idea of possible photochemical reaction 
mechanisms and reaction branching ratios (j52y i9 1 [53 t[5H) . 

4.2.3 The Challenge of Excitonics To judge by the numerous articles 
which have appeared in the Annual Review of Physical Chemistry, excitonics, 
or the study of excitons, is of great and recurring interest to our community. 
It would most likely be described differently by a solid-state physicist and by a 
physical chemist or biophysicist. A solid-state physicist would probably begin 
with the image of delocalized states with the expected broad absorption spectra. 
Sharp spectral features are then often associated with localized excitations which 
could propagate throughout a material and eventually the electrons and holes 
could separate to do something interesting such as generate electricity or drive 
a chemical reaction. A physical chemist or biophysicist would probably begin 
by dissecting an extended system, such as a molecular solid or a biomolecule, 
into components which are likely to have localized excitations and then use an 
excitonic model to try to understand how excitations interact to lead to spectral 
shifts or energy transfer. In the end the concepts are not very different except 
that the solid-state physicist is localizing a delocalized phenomenon while the 
chemist is delocalizing localized phenomena. 

These two points of view are now frequently coming together in the area of 
nanoscience where solids and molecules frequently meet at the nanointerface. 
This is particularly apparent when modeling dye-sensitized solar cells (j55p . There 
is a level matching problem in these systems where the excited-states of dye 
molecules must be aligned with the band structure of transparent semiconductor 
particles so, for example, that electrons can pass into the conduction band when 
neutral excitations break up at the molecule-crystal interface. Calculations such 
as that of Ref. (j56p give an idea of what is currently possible when modeling 
dye-sensitized solar cells with TD-DFT. 

Excitons are also highly relevant in photobiology where exitations in nearby 
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proteins are coupled, thus allowing energy to be captured by antenna systems 
and transported to active sites where electronic energy is used to do biologically- 
important chemistry. Both for practical and conceptual reasons, it makes sense 
to use subsystem density-functional theory. Johannes Neugebauer has been par- 
ticularly successful in this regard in implementing and improving the subsystem 
version of TD-DFT ([SZllSl]) initially proposed by Casida and Wesolowski (f27|) . 
A recent application of this theory is to the light harvesting complex II of green 
plants (1591) . 

Given the system sizes involved and the importance of also simulating exciton 
dynamics, it is probable that the TD-DFT will have to be combined with semi- 
empirical models, such as time-dependent DFTB (|60p. at some point in order to 
describe the full complexity of excitonic phenomena. 

4.3 Problems, Problem Detection, and Solutions 

By discussing the applications of TD-DFT before discussing the problems of TD- 
DFT, we may have given the false impression that classical TD-DFT is a blackbox 
method. It is far from that! Just as pure Kohn-Sham DFT looks like the HF 
(or rather the Hartree) approximation but behaves very differently, TD-DFT 
even if put in the form of the classic RPA equation of quantum chemistry does 
not behave like TD-HF. Conventional LR-TD-DFT in its basic AA using pure 
DFAs often, but not always works better than LR-TD-HF. It has been found to 
work best for (i) low-energy (ii) one-electron excitations involving (iii) little or 
no charge transfer and (iv) which are not too delocalized. This section examines 
the reasons for these limitations and gives a few suggestions as to what can be 
done when a better solution is needed. Interestingly many of the problems are 
not due to the adiabatic approximation so that our analysis may be carried out 
in terms of Jacob's jungle gym (Fig. [3]). Discussion of problems associated with 
the adiabatic approximation and possible solutions will be delayed till sections [5] 
andm 

4.3.1 Problems associated with E^c There is an intimate relationship 
between the ground- and excited-electronic states in conventional TD-DFT. Not 
only are the excited-state PESs obtained by adding TD-DFT excitation energies 
to the DFT ground-state energy, but the same A and B matrices that occur in 
LR-TD-DFT also occur in a stability analysis of the DFT ground state (161( 113 ^ 
[22]) . In fact the theory closely parallels what happens in HF stabililty analysis 
(SecESI). 

In order to understand the similarities and differences between what happens 
in DFT and HF, let us first look at the classic example of H2. In Hartree- Fock 
theory, stretching the bond beyond the Coulson-Fischer point leads to symmetry 
breaking in the sense that a DODS solution becomes lower in energy than a SODS. 
The reason is well understood in wave- function theory and has to do with presence 
of ionic contributions in the DODS HF wave function. The situation is different 
in DFT. Assuming NVR, we must have a SODS solution. There is thus a unique 
nodeless orbital, il^{r), which is found by taking the square root of the density. 
The noninteracting potential is then calculated as, fs(r) = (V'^ip{r)/2'ip{r)) — e, 
with the orbital energy, e, adjusted so that Vg goes to zero at infinity. The NVR 
assumption may then be explicitly verified by solving the Kohn-Sham equation 
to make sure that e does indeed correspond to the lowest orbital energy of the 
noninteracting system. In this way, it is seen that symmetry breaking does not 
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occur in the exact theory. 

But it does occur for DFAs even though it is not as severe as in HF. In particular 
the Coulson-Fischer point for pure density functionals occurs at larger H2 bond 
distances than in HF, but it does occur and is recognizable because the square 
of the lowest triplet excitation energy, ujj^, goes to zero exactly at the Coulson- 
Fischer point before becoming negative (an imaginary triplet excitation energy) 
at larger distances (jlSh . Incidently it is perfectly possible for an SCF calculation 
beginning with a SODS guess to converge to a SODS solution even when there 
is a lower-energy DODS solution ()62p . In this sense, a TD-DFT calculation may 
serve as an important check on the correctness of the ground-state calculation. 
However when the goal is to calculate excited-state PESs, the presence of triplet 
and near singlet instabilities is highly inconvenient. As in the HF partial 
solution is then provided by the TDA whose excitation energies "inherit" some 
of the variational nature of CIS and so do not collapse when the ground-state 
solution becomes unstable p3y i8 1[T9|) . 

A second phenomenon which can occur when breaking a bond is the occurance 
of a cusp in the ground-state PES when orbital fillings suddenly change ps y GS p . 
The wave-function solution to this problem is to allow configurations to mix and 
so have the wave function be a linear combination of determinants corresponding 
to the two different orbital fillings. In contrast this usually corresponds in DFT 
calculations to the situation of effective failure of NVR when the LUMO falls 
lower in energy than the HOMO. Since most DFT programs enforce Aufbau 
filling of the orbitals, the LUMO is immediately filled so that it becomes the next 
HOMO and the old HOMO the new LUMO. However the LUMO remains lower 
than the HOMO and the calculation simply fails to converge (llSpiQp . As already 
discussed in Sec. 13. 1^ this "strange" orbital filling means that we are, in fact, in 
the presence of a spin-wave instability. The conventional DFT solution is to go 
to ensemble theory and allow the orbitals to have fractional occupation number. 
Unfortunately it is really far from clear how well conventional DFAs apply to the 

ensemble situation. 

This situation is probably one of the contributing factors for the development of RDMFT: Reduced density- 
reduced density-matrix functional theory (RDMFT). In principle, the density is matrix functional theory. 

replaced with the 1-RDM. However in practice it is the natural orbitals and their 

corresponding occupation numbers which become the fundamental variables. As 
such it is very much like ensemble density functional theory but also has the 
possibility of handling nonlocal potentials (i.e., one-electron operators which are 
not simply multiplicative functions). Exact equations may also be derived for 
two-electron systems since Lowdin and Shull showed that the 2-electron wave 
function may be expanded exactly (jMp . The spatial part is given by. 



where the natural orbitals, (pi{r), are assumed to be real, their complex phase, 
Oj, has been exhibited explicitly, and the associated natural orbital occupation 
numbers are denoted rij. Defining the phase is the major problem in the two- 
electron theory, but typically the minority natural orbitals are found empirically 
to have a negative phase with respect to the dominant phase. Very good exact 
functionals may then be developed for the two-electron case which, for example, 
give the correct dissociation of H2 . Notable early work attempting to develop 1- 
RDM functionals for systems with more than two electrons include that of Miiller 




(38) 
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()65p and Goedecker and Umrigar (I66p . For more recent work see, for example, 
Ref. (j67p and references therein. 

Finally there is a problem with conical intersections which arises precisely 
because excited-state PESs are calculated in TD-DFT as, Ej = Eq + ooj. For a 
molecule with Nf internal degrees of freedom, Q, the PES is an A^^'-dimensional 
object in an [Nj + l)-dimensional space (because we have added the energy 
axis). When two PESs cross, we have the constraints (i) i?/(Q) = Ej{Ql) and 
(ii) Hjj{Q) = 0. Constraint (i) is always present in TD-DFT and reduces the 
dimensionality of the crossing space to a hyper line always. Constraint (ii) should 
further reduce the dimensionality of the intersection space to a hyperpoint — i.e., 
a conical intersection where there are two independent directions going away from 
the intersection space in hyperspace where the surfaces separate (unless of course 
we have a diatomic so that A'^^ = 1 in which case there is an avoided crossing 
instead.) However Levine, Ko, Queenville, and Martinez have pointed out (j68p 
that constraint (ii) is absent for interactions between the excited states and the 
ground state because of the use of the formula Ej = Eq + ojj. We are thus left 
with the condition that conical intersections cannot exist in TD-DFT. While this 
is certainly correct, in practice they may be approached closely enough to make 
mixed TD-DFT/classical surface hopping photodynamics calculations at least 
qualitatively useful for investigating photochemical reaction mechanisms (jl9p . 

4.3.2 Problems ASSOCIATED WITH The problems arising from E'l^c typ- 

ically arise at molecular geometries corresponding to bond making or breaking in 
the ground state. Problems associated with Vxc occur even in the Franck-Condon 
region near the ground-state equilibrium geometry. In Sec. 14. 1^ we argued that 
the Kohn-Sham orbital energy difference should lie between the corresponding 
singlet and triplet excitation energies [Eq. (I37p ] and that all three quantities will 
be well approximated by the Kohn-Sham orbital energy difference for Rydberg- 
type excitations. Since the limit of a given Rydberg series is just the ionization 
potential of an orbital tpi, we see immediately that the TD-DFT ionization thresh- 
old is at minus the HOMO energy. However practical DFAs have minus HOMO 
energies which are typically several electron volts too low compared to the exper- 
imental ionization potential This means that finite basis TD-DFT calculations 
lead to a collapse of the higher excited states in the region between minus the 
HOMO energy and the true ionization potential (j69p . The reader needs to be 
keenly aware of this point because the collapse may not be obvious when only 
medium-sized basis sets are used, though it will occur once large-enough basis 
sets are employed. 

One approach that could be taken is simply to ignore the problem. After all, 
one might think that that the oscillator strength distribution should be approxi- 
mately correct even above the TDLDA ionization threshold at —enoMO (|70p . Of 
course, transitions to bound states will no longer be sharp, making assignments 
more difficult, but the spectral features will often still be present. 

Probably a better approach is to correct the asymptotic behavior of Vxc by 
a shift-and-splice approach (see Ref. (j7ip and references therein.) This is suffi- 
cient, for example, to obtain qualitatively reasonable Rydberg PESs and hence to 
see valence-Rydberg avoided crossings in TD-DFT calculations for formaldehyde 
(H2C=0) (J72D. An even more sophisticated approach is the previously-mentioned 
OEP model. 

A subtler issue which is only occasionally discussed in the DFT literature is 
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that the relative orbital energies of occupied orbital energies obtained with DFAs 
may not be correct compared with the orbital energies which should come out 
of exact Kohn-Sham DFT calculations. In Refs. (j22t[73]) OEP calculations were 
used to show that errors in ethylene (H2C=CH2) asymptotically-corrected LDA 
TD-DFT excitation energies were due to small relative errors in the occupied a 

and vr MO energies. 

A final issue which is far from minor is the problem of charge transfer excita- RSH: Range-separated 

tions (|74p . Up to this point we have argued that DFT orbitals are preprepared hybrid. 

to describe excitations, but a better statement is that they are preprepared to 

describe neutral excitations. Indeed, as discussed just after Eq. (j20p HF theory 

seems to be very well designed for describing long-range charge transfer. An 

inspection of the TDA-TD-TDDFT formulae in Table [1] shows that the long- 
range charge transfer energy in TD-DFT is expected to be just the orbital energy 

difference, €a — €i, which is badly underestimated because of problems with the 

asymptotic behavior of Vxc related to the PNDD. This also suggests that the 

overlap between the initial and final orbitals in an excitation could be used as a 

diagnostic tool to anticipate when TD-DFT will break-down due to charge trans- 
fer excitations ()75p . As the linearized ASCF DFT formulae in Table [1] show, 

this problem is not as severe in the ASCF approach which takes partial account 

of orbital relaxation energies and this has been the basis of one proposal for a 

charge transfer correction psp . That paper also makes the point that not every 

charge transfer excitation is problematic in TD-DFT — only the ones involving 

density relaxation. In fact, MuUiken's prototypical charge transfer excitation 

in H2 is an excellent example where density relaxtion is small and TD-DFT has 

little problem. However a better answer can be found in the range-separated 

hybrid (RSH) approach (see e.g., Ref. ([76ll and references therein.) The basic 

idea here is to separate the electron repulsion, l/ri2, into a short-ranged part, 

(l/ri2)g^ = erfc(7ri2)/ri2, which is described by a pure DFA and a long-range 

part, (l/ri2)^^ = erf(7ri2)/ri2, which is described by a suitable wave-function 

theory such as Hartree-Fock. There are several variations on the name and precise 

functional forms used in this theory, but they do yield a dramatic reduction in 

errors due to the underestimation of charge transfer excitations. Drawbacks are 

the presence of a possibly system-dependent parameter (7 in the above example) 
Interestingly a recent OEP study shows that errors in the charge transfer ex- 
citation of HeH^ are not due to errors in the xc-kernel, but rather to errors in 

Kohn-Sham orbital energies arising from DFAs (|77p . A second OEP study con- 
firms that a significant part of the errors in the adiabatic TD-DFT is due to the 

use of DFAs, but goes on to show that an important residual error is due to neglet 

of the frequency dependence of the xc-kernel ([Tsp . 

4.3.3 Problems associated with f^c There are two main problems asso- 
ciated with fxc- The first has to do with the need to include not just Iplh exci- 
tations but also higher excitations in order to treat some problems. For example 
explicit 2p2h excitations are needed to describe the first singlet excited state of 
butadiene, CH2=CHCH=CH2, has significant double excitation character ()79p . 
Also explicit higher-electron excitations are needed for a proper description of 
excitations of molecules with open-shell ground states ()80p . Although exact TD- 
DFT should describe this correctly, making the adiabatic approximation restricts 
the number of solutions of the LR-TD-DFT equation to the dimensionality of the 
matrix, namely only the number of Iplh excitations and de-excitations. Including 
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some frequency-dependence in f^c could allow a nonlinear feedback mechanism 
allowing the matrix equation to have additional solutions. At one point, it was 
hoped that going beyond LR to look at, say, simultaneous two-photon absorption, 
would allow two-electron excitations to be treated within the adiabatic approx- 
imation ()8ip . however it is now clear that this is not the case. In particular, 
the poles of the dynamic second hyperpolarizability are identical to the poles of 
the dynamic polar izability (182p . which is to say the one-electron excitations of 
adiabatic LR-TD-DFT. Interestingly two-electron excitations are accessible via 
sequential absorption in real-time TDHF and TD-DFT ()83p . 

There is a second problem associated with fxc, but also with Vxc and even 
Exc- This problem goes by different names in the literature. We will refer to 
it as the "scale-up catastrophe," because it arises when a system is made larger 
and larger. Now, an important goal of DFT is to extrapolate ab initio accuracy 
to molecules too large to otherwise treat by conventional (i.e., HF-based) ab 
initio wave- function calculations. The usual recommendation is to test DFT 
calculations on small molecules where comparisons can be made with accurate ab 
initio calculations and then to assume that DFT calculations which are reliable 
for a given class of molecule and molecular properties will remain accurate as 
the size of the molecules treated is increased (i.e., scaled-up). When instead 
inaccuracies increase as the size of the system is increased, then we have a scale- 
up catastrophe. A perhaps little-known example of a scale-up catastrophe in LR- 
TD-DFT occurs in periodic systems where the xc-kernel vanishes on the lower 
rungs of Jacob's ladder due to improper scaling with respect to the number of 
K points. ()84p The result is that the first vertical excitation energy reduces to 
just the difference between the HOMO and LUMO orbital energies, which is 
too low. Nevertheless the spectrum may remain essentially correct because the 
oscillator strength of this transition may be negligeable. The better known and 
very disturbing example of a scale-up catastrophe is that studies on oligomers of 
conducting polymers showed that dynamic polarizabilities could be overestimated 
by arbitrarily large amounts by simply going to a large enough oligomer (j85ll86p . 
(I85p 87 t[86]l There are different ways to try to understand this error. For example, 
most DFAs are based upon the homogeneous electron gas which is a model for a 
perfect conductor. Applying an electric field to a conductor leads to the build up 
of surface charges which cancel the applied fields in the interior of the conductor 
(Faraday effect). If the system is behaving too much like a conductor, then the 
HOMO-LUMO gap is too small due, for example, to an improper description of 
the PNDD. That is an explanation mainly in terms of orbital energies and hence 
Vxc- However the reaction field which cancels the applied field in the interior of 
the conductor is given by the response of the self-consistent field which depends 
on the description of the xc-kernel. This xc-kernel must be "ultra-nonlocal" in 
the sense that it must react to charges built up on the surface of the conductor. 
On the lower rungs of Jacob's ladder, fxc'"^ (^i, r2), is roughly diagonal in (ri, r2). 
This is in contradiction with a rough estimate, exact in the case of two-electron 
systems, of the exchange-only part of the xc-kernel which suggests that it be 
more like — |7(ri,r2)p/ (/9(ri)ri2/o(r2)) ([88]) . This means that a diagonal fxc is 
certainly incorrect. There are few choices for dealing with this problem — either 
one restricts oneself to medium-sized systems before the onset of the scale-up 
catastrophe or one climbs Jacob's ladder to use more sophisticated functionals. 
One way to restrict oneself to medium-sized systems is to use a subsystem theory 
such as the one used to treat excitons in light harvesting complex II, so that TD- 
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DFT is only restricted to medium-sized molecules at any given moment. Climbing 
Jacob's ladder to functionals with contributions from Hartree-Fock exchange is 
also helpful as is the use of RSHs (|89p90p9H92p . Another approach which is helpful 
in some cases is time-dependent current-density functional with the Vignale-Kohn 
functional (193p45p . This is because the current can carry information about 
perturbations at a distance. 



5 GOING BEYOND THE ADIABATIC APPROXIMATION 

The exact TD-DFT potential [Eq. (j33p ] is a functional of the whole history of 
past densities and the initial conditions. These memory effects, are completely ne- 
glected in common TD-DFT applications using the AA. In fact, the AA assumes 
that the xc-potential reacts instantaneously and without memory and so may be 
expressed in terms of the ordinary DFT xc-energy evaluated for the density at 
that time, 

Thus the AA xc-potential is only a functional of the function p^„{t) of the three 
independent variables r = (x, y, z) at fixed time, and not a functional of the 
full function p„{r,t) of four independent variables, a fact that we have tried to 
emphasize by our choice of notation. In this sense, TD-DFAs using the AA 
work fairly well when the non-interacting system is a reasonable physical first 
approximation to the interacting system, but fail dramatically in cases like 2p2h 
excited states ([MI), conical intersections ([68l[T8]) . charge transfers and Rydberg 
excitations (fTillTSj) . excited states along the bond dissociation coordinate (|95|) . 
one-dimensional extended systems (I96 p97p . band-gaps in solids (|98ll and spectra 
of semiconductors ([98l [99 |[T00]) . In such cases, the contribution of past densities 
to the xc energy becomes essential, requiring memory functionals specifically 
designed for TD-DFT. In this section, we will review the most successful memory 
functionals available at the present time. Few attempts have tried to approximate 
directly the xc-action (llOip . Instead, most memory functionals approximate the 
derivatives. Therefore, our discussion will focus on the construction of memory 
xc-potentials and specially memory xc-kernels. 

Memory introduces strong requirements for TD-DFT functionals ()102p . which 
forbids the intuitively simple functional approximations of ground-state DFT. 
For example, a simple parameterization of a spatially-local memory LDA xc- 
kernel for the homogeneous electron gas (jlOSp violates, by construction, the zero- 
force theorem (13Hll04p . introducing unphysical damping effects during the time 
propagation of the density and a corresponding lose of total energy. Several 
schemes have appeared that enforce the sum rules a posteriori (jlOSyiOG"]) . Most 
importantly, these studies showed that memory functionals must be non-local in 
space. This property forbids the application of the popular gradient expansion 
to memory functionals (jl07p . ultimately a consequence of the ultra- nonlocality 
property of the exact xc-functional. 



5.1 Current and Lagrangian Density- Functional Theory 

The ultra-nonlocality property of the exact functional indicates that the density is 
perhaps not the best reduced variable to express time-dependent xc-effects. This 
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redirected the attention to other time-dependent theories where memory effects 
are described by local functional, such as in TD-CDFT pU7lllU8l[T09l[TT0l[T0T]) . 
in Lagrangian TD-DFT (|1111I1121[TT3]1 . 

TD-CDFT is very much related to TD-DFT, though its range of applications 
is broader. In fact, the one-to-one correspondence between currents and exter- 
nal potentials is the first step in the Runge-Gross existence proof for TD-DFT 
([3]). However, the TD-CDFT proof does not resort to any surface integral, which 
makes its application to solids more straightforward than TD-DFT. The first 
memory functional for TD-CDFT is due to Vignale and Kohn, which consists 
of a parameterization of a memory LDA xc-kernel tensor for the homogeneous 
electron gas ()107p . This kernel is local in space and, unlike the corresponding 
in TD-DFT, it satisfies the zero-force theorem. However, the kernel introduces 
too strong damping effects that induce a decoherence mechanism, which allows 
the density to decay back to a ground-state with higher entropy, different from 
the real ground-state, in which the total energy is conserved pi4p . The Vignale- 
Kohn kernel has been used successfully to calculate the static polarizabilities of 
one-dimensional conjugated chains, which were largely overestimated by common 
adiabatic LDA ()96p . However, it fails for one-dimensional hydrogen chains, re- 
quiring more sophisticated kernels ()97p . The Vignale-Kohn kernel also shows the 
characteristic double-peak spectrum of bulk silicon (jllSp . a feature that is not 
present in the adiabatic LDA spectrum. Memory kernels of TD-CDFT can then 
be map into TD-DFT kernels maintaining part of the features of the Vignale- 
Kohn kernel (jllOp . 

Lagrangian TD-DFT has less applications to realistic systems, but it has pro- 
vided deeper insight on the role of memory. Lagrangian TD-DFT is a reformu- 
lation of TD-DFT, in which the electronic density is re-expressed in a reference 
frame (Lagrangian frame) that moves along with the fluid (Illipil2]) . This allows 
the memory to be expressed locally in terms of a the position of fluid elements 
and the deformation tensor. The deformation tensor accounts for the Coulomb 
coupling of a differential volume of fluid at position r and time t with was at a 
different position r' at an earlier time t' < t. Both TD-CDFT and Lagrangian 
TD-DFT are equivalent in the linear regime (I113p . 

5.2 Optimized Effective Potential Approaches 

All the aforementioned memory functionals lack a systematic route to improve 
upon the present approximations. A more systematic way to construct mem- 

ory functionals is to use the TD-OEP, in which many-body perturbation theory 

MBPT: Many-body per- (MBPT) quantities are mapped into TD-DFT (jll7[lll8[[Tl^ [T2n l[T2Tl[T22l[T23ll . 
turbation theory. An easy way to establish this mapping is to make use of the Kohn-Sham as- 

sumption, i.e., requiring that the interacting time-dependent density equals the 

non-interacting one (jll9p . This is the basis of the time-dependent Sham-Schliiter 
equation that allows us to relate the self-energy to the xc-potential and the kernel 
of the Bethe-Salpeter equation to the xc-kernel. In this way, consistent approx- 
imations are possible, that guarantee the satisfaction of exact conditions (|12ip . 
There are many ways to derive the time-dependent Sham-Schliiter equation, but 
similar functionals are obtained if the same MBPT approximation is used. For 
deriving xc-kernels, a simple way is to require that the diagonal of the MBPT 
response functions leads to the same spectrum as the TD-DFT response function 
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[Eq. ([22])], which leads to 

fUr,r';u)= (40) 
J d^n J (fr2 J (frs J (i^r4A(r;ri,r2;a;)i^:(ri,r2;r3,r4;a;)Aj(r3,r4;r';u;) , 

relating the exact TD-DFT xc-kernel to the kernel K{uj) = IlJ^{io) — n~^(a;) 
defined from the polarization propagator from MBPT (|123p . The mapping is 
possible thanks to the localization function A, which is expressed in terms of the 
response functions [Eq. ([22]) ]. 

A(ri;r2,r3;a;) = x'^ri, r2; w)x(r2, rs; (j) . (41) 

The localizer A^, appearing in Eq. (|40|) . is defined as in Eq. (|41|) except in terms of 
the noninteracting response function (^^(a;)). The localizer introduces an extra 
frequency-dependence in the kernel, which brings in the MBPT quantity the 
ultra-nonlocality property of the exact TD-DFT functional (|124lll22l[T25l[T23]l . 
To date, little is known about the exact role of the localization step, though it is 
known to have interesting properties. For example, Gonze and Scheffler showed 
that the single-pole approximation when one applies the linear Sham-Schliiter 
equation in Eq. (jlT]) gives ()124p 

ujia = Aeia + {ia\f^c{Aeia)\ai) = Aej^ + {ii\K{Aeia)\aa) , (42) 

where Aeia = £a — £«• In fact, a cancellation of the localizer is observed, and the 
exact xc-kernel is equal to a different matrix element of the kernel derived from 
MBPT quantities. A general treatment of localization is plagued with numer- 
ical problems, due to the difficulties in inverting the singular density response 
functions that are involved (jl26yi27"j) . In the linear-response regime, an interest- 
ing work around can be used to fully account for the localization effects without 
the need of explicitly inverting any response function. This is achieved by re- 
expressing LR-TD-DFT as the response of the non- interacting potential (jl28p . 

Using a first-order approximation of K, one obtains the TD-EXX kernel (|1181 
I120yi21pl23p . This kernel provides correct description of charge transfer (|78p and 
correct position of the continuum states ([6|) due to the correct asymptotic decay of 
exchange. Although the exchange kernel in MBPT is frequency independent, the 
full TD-EXX is a frequency-dependent functional (jll8lll27P due to the localizer 
A. Therefore some (minor) memory effects are taken into account ()129p . In 
general, the frequency dependence of the localizer in the TD-EXX is not enough 
to include 2p2h- and higher-ph excited poles (jl30U13ip . and frequency-dependent 
MBPT kernels are required pJ lfTH2|fT2HlfT33l) . 

The lack of 2p2h excitations is an endemic problem of the adiabatic approx- 
imation dHj). These poles should be present, since TD-DFT is formally exact. 
However, they are not present when an adiabatic xc-kernel is used in the LR- 
TD-DFT equations. Clearly, the number of solutions of Eq. ([36]) is given by the 
dimension of the matrices, that is, the number of Iplh states. Recovering of 
2p2h states require a frequency-dependent xc-kernels (1134p94p . which make the 
LR-TD-DFT equations non-linear and extra solutions appear. Frequency depen- 
dence can be included via the solution of a Scham-Schliiter equation. The MBPT 
kernel of the Bethe-Salpeter equation has to be constructed carefully, since oth- 
erwise the resulting xc-kernel might introduce spurious transitions (jl32pi33"]) . 
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The first successful tj-dependent xc-kernel is the kernel of dressed TD-DFT, 
first proposed by Burke, Maitra and coworkers (|1341[M|) to include one 2p2h at 
a time, and later generalized to any number of 2p2h states by using full-featured 
MBPT techniques (j80pi23p . The kernel of dressed TD-DFT can be viewed as a 
hybrid kernel between the adiabatic xc-kernel and the frequency-dependent part 
of the MBPT kernel, 

fUr, r'; = /,1.^(r, r') + i^-(r, r; r', r'; oj) , (43) 

where K^(r,r;r' ,r';Lj) corresponds to the frequency-dependent part of the ker- 
nel derived from a second-order polarization propagator approach (for the math- 
ematical expression see ()123p and ([6]).) Several variations on dressed TD-DFT 
are possible. However Ref. ([6]) explores quite a few of these and comparison of 
the performance of dressed TD-DFT against benchmark ab initio results for 28 
organic chromophores goes far towards establishing the best protocol to use. 

5.3 Time-Dependent Reduced Density Matrix Functional The- 
ory 

Important recent work has been done by Giesbertz and coworkers extending 
RDMFT to the time-dependent case ((133 [95l [Ml [EZl [l38l [IMl EQ]) . Oddly 
enough most of these attempts fail to obtain the correct w — ?• limit of the 
adiabatic approximation for TD- RDMFT even for two-electron systems ()139p . 
However including not only the orbitals and occupation numbers in the time- 
dependent response equations but also their phases allows exact solutions to be 
obtained for two-electron systems ()138p . The observation of the importance of 
phases is reminiscent of our observations regarding spin-wave instabilities and 
violation of NVR and might be taken as further proof of the need to include 
phase information TD-DFT. Note that this phase information is not included in 
the 1-RDM, so that a new functional dependence is needed which goes beyond 
the density-matrix alone. It may perhaps be assimilated with a memory effect. 

6 GOING AROUND THE ADIABATIC APPROXIMATION 

In the last section we described progress in TD-DFT going beyond the adiabatic 
approximation but trying to keep the essential spirit of TD-DFT. In this section 
we describe approaches which try to solve the problems of the TD-DFT adiabatic 
approximation while still retaining something resembling frequency-independent 
kernel. They may be thought of as ways to go around the problem. 

6.1 Reconciling TD-DFT and ASCF 

We have seen that the ASCF method and its MSM variant predate TD-DFT 
in quantum chemistry. The older method is particularly well justified for cal- 
culating the first ionization potentials and electron affinities, the lowest triplet 
excitation, and perhaps also transitions in general to states dominated by a sin- 
gle configuration. Indeed the numbers produced by the MSM ASCF method and 
TD-DFT are often quite similar, though the formulae are different. It has been 
a puzzle since the beginning of applications of TD-DFT in quantum chemistry 
dU to try to understand why this is so. Reconciling the two methods might also 
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allow some of the better aspects of the MSM ASCF method to be used in TD- 
DFT, namely the ease of handling 2p2h and higher-ph excitations in the ASCF 
method or apparently better ability to describe density-relaxation during excita- 
tion. A first step in this direction were presented by one of us in 1999 and may 
be found in the on-line proceedings (jl4ip . These ideas were later used to develop 
a charge-transfer correction for TD-DFT (jl3p . Hu, Sugino, and Miyamoto have 
followed up this work to develop their own approach to correcting problems with 

underestimated Rydberg and charge transfer excitations. (jl42yi43"]) 

More recently these ideas have been very nicely followed up in Tom Ziegler's CV: constrained varia- 
group in Alberta in the form of constrained variational DFT (CV-DFT) ()144pi45^ tional 

I146pi47|) . The energy is to be minimized with respect to the transformation, (/)^ = 

'Pi + J2a Ua,i4>a + 0^'^\u), of the Kohn-Sham orbitals, subject to the constraint 
WU = 1. The second-order CV(2)-DFT energy expression is already familiar 
from stability analysis (Eq. (I2ip with A = 1). Minimizing to this order just gives 
the TD-DFT TDA. However fourth order CV(4)-DFT with relaxed orbitals gives 
a much improved treatment of charge transfer excitations with out inclusion of 
any nonadiabatic frequency dependence (|148p . 



6.2 Spin-Flip TD-DFT 

Another way to go around the adiabatic approximation — this time to access a 
2p2h excitation — is to use the spin-flip method. The basic idea is illustrated in the 
TOTEM in Fig. [H To do this properly, we must consider the Kohn-Sham orbitals 
as spinors, (■0Q(r), '(/'^(r)), which can rotate their spin-direction at different points 
in physical space. The density also acquires spin-labels. 



P(r) 



(44) 



So the xc-kernel, /^^''^^''^^'^^(r, r'), also acquires additional labels. However this 
is just for the purposes of deriving a more general LR-TD-DFT equation. In the 
end, we will evaluate the xc-kernel using the usual collinear Kohn-Sham theory 
(i.e., each spin orbital is either a or /? but not some linear combination of the 
two). However it turns out that conventional pure DFAs do not allow spin- flips. 

In order to get good agreement with experiment. Ana Krylov and coworkers 
(jl491ll50j) found it necessary to use a significantly higher amount of HE ex- 
change (50%) than is typically used for ground state properties (~ 25%). Even 
higher percentages of HE exchange (> 50%) have been reported to be necessary 
for calculating second hyperpolarizabilities of diradical systems by this spin-flip 
method ()15ip . Although the use of a different functional for ground and excited 
states is disturbing, the basic idea is admirable and this method continues to 
be used (jl5Hll52]ll53p . In particular, this is the spin-flip TD-DFT approach 
mentioned in the introduction in the context of its recent use by Minezawa and 
Gordon who found the method to give a relatively good description of conical 
intersections in ethylene ()154p . 

The next and most recent major advance in spin- flip TD-DFT came with an 
article by Wang and Ziegler ()155p . It is intimately related to work by Wenjian Liu 
and coworkers on relativistic four-component TD-DFT ()156p . Wang and Ziegler 
proposed that any pure spin-density xc-functional, Exc[pa, P/3], could be used to 
make a noncollinear xc-functional suitable for spin-flip calculations by making 



30 



Casida k, Huix-RotUant 



the substitution, 



Pa^ P+ = -{p + s) , pp^ p^ 



(45) 



involving two quantities which are invariant under a unitary transformation of 
the spin coordinates. These are the total charge density, p = pa,a + Pp,i3, and the 
magnetization, s, whose square is given by, = {pa,a — P(S,I3)'^ + 2 (^p'^ /? + /'I «) • 
The collinear limit of s is just the spin-polarization, s — t- Pa — Pi3, after an 
appropriate choice of phase. The factor of 1/2 has been introduced by us ([63|) so 
that, 

P+^ Pa , p-^ Pfi, (46) 

in the same limit. After taking derivatives and the noncollinear limit, the xc- 
kernel becomes, 



J xc J xc 

J xc 

J xc J xc 



xc J xc J xc J xc 

•I3,l3;a,a fli,fi;P,l3 fli,fi;a,fi fl3,^;[S,a 

xc J xc J xc J xc 

xc J xc J xc J xc 

•0,a;a,a ffi^a^fif} fl3,a;a,l3 fl3,a;l3,a 

xc J xc J xc J xc 



J xc J xc 



I xc 






' xc 











„/3 



Pa-P/3 












(47) 

This approach to spin-flip TD-DFT has been applied to the dissociation of 
H2 ()155p and to calculate the spectra of open-shell molecules (I157pi58lll59p . 
We have applied the Wang-Ziegler approach to the photochemical ring open- 
ing of oxirane (H2COCH2) ()63p . While this approach definitely overcame the 
cusp problem for C2V ring opening and does allow coupling of the ground- and 
excited-states so that there is a true conical intersection, the position of the con- 
ical intersection was found to be intermediate between the CIS seam and the 
complete active space SCF conical intersection (j63p . Very recently analytical 
derivatives have been worked out for the Wang-Ziegler method ()160p . 

The Wang-Ziegler approach has also been proposed as the basis of a more 
general spin-coupled TD-DFT (116ip and applied (|162p . Very recently it has 
been used to solve the problem of spin-adaption for LR-TD-DFT in open-shell 
molecules (|163pi64]) that had been highlighted by earlier work on spin contamina- 
tion in TD-DFT (|165p . Hybrid spin- flip functionals are also beginning to appear 
(1162111661) • 



7 PERSPECTIVES 



In coming to the final section of this review, we are deeply conscious of the many 
aspects of TD-DFT which we have neglected and the shallowness of our treatment 
of those aspects of TD-DFT which we have been able to treat in this small 
overview of what has become a rather vast field. Nevertheless we hope that the 
reader has appreciated the historical interplay between DFT, TD-DFT, and the 
development of DFAs. It is a fact that adiabatic TD-DFT in seeking to calculate 
new types of properties, has placed brought into sharper focus known problems 
of and has placed new demands on DFAs. This in turn has pushed further 
improvement of DFAs for classical DFT and pushed more detailed exploration of 
Jacob's ladder. 

Also significant has been the difficulty of answering the Runge-Gross challenge 
to go beyond the AA and include memory in TD-DFT. Section [5] has discussed 
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methods for trying to go beyond the AA while still trying to preserve the true 
spirit of (TD-)DFT — that is, to do as much work analytically before asking the 
computer to do the work. In the end, this means taking things from MBPT and 
putting them into the functionals. When we do the opposite — that is adding 
MBPT on top of the functionals — then that is an indication of something that 
we do not yet understand well enough and need to work on some more. 

The penultimate section (Sec. [6]) discussed primarily pragmatic ways to try to 
circumvent problems created by the AA without abandoning the AA. These, plus 
dressed TD-DFT (Sec. [5|), are the tools we have now but which may merge in the 
future with the more rigorous approaches of Sec.[5]to give something approaching 
an answer to the Runge-Gross challenge. 

Only time will tell . . . but then time is such an important and fundamental 
physical parameter! 

Summary Points 

1. After almost three decades, the Runge-Gross theorem is largely uncontested 
as the foundations of a rigorous theory. Initial problems in the definition 
of the action have been heavily invetigated, leading to increasingly satisfac- 
tory definitions of appropriate xc-actions. Now most fundamental TD-DFT 
research is devoted to improving TD-DFAs. 

2. TD-DFT is making good progress towards a general application to photo- 
chemical problems. Charge transfer, Rydberg excited states and 2-particle/2- 
hole states can now be correctly treated within LR-TD-DFT. The treat- 
ment of conical intersections will require overcoming the non-interacting 
tJ-represent ability problem in DFT, which requires fractional occupation 
numbers and very possibly complex orbitals. 

3. In general, LR-TD-DFT using the Tamm-Dancoff approximation leads to 
better excitation energies and it should be preferred when calculating po- 
tential energy surfaces. However, the Tamm-Dancoff approximation can 
lead to inaccurate oscillator strength distributions. 

4. Most of the problems attributed to LR-TD-DFT can be simply resolved 
by using decent approximations of the DFT functionals. From practical 
experience, asymptotically-corrected hybrid functionals with 20-25% of HF 
exchange lead to the best results for singly-excited states of finite systems. 

5. Memory functionals have extended the range of applicability of TD-DFT, 
though some problems are yet to be resolved. The longstanding problem 
of 2-particle/2-hole states can be resolved by adding to the adiabatic ap- 
proximation a frequency-dependent component derived with the help of 
MBPT techniques. Though the mixture requires careful application to 
avoid double-counting of correlation, present applications show encourag- 
ing results. 

6. Spin-Flip LR-TD-DFT currently offers a practical solution for photochemi- 
cal problems. It adds some double excitation character and it avoids regions 
of the potential energy surface with noninteracting u-representability prob- 
lems, offering an easy route to calculate conical intersections. Also, it can 
successfully treat excited states of some open-shell systems. 



32 



Casida k, Huix-RotUant 



Future Issues 

1. Even though some advances have been made, memory functionals are still in 
their infancy, and better functionals are necessary for extending the range 
of application of TD-DFT. It seems more appropriate to think in terms 
of a different Jacob's ladder for memory functionals in TD-DFT with the 
rungs defined according to the functional dependence as: TD-DFT (den- 
sity), TD-CDFT (density and current) and L-TDDFT (fluid position and 
deformation tensor), TD-OEP (orbitals), TD-RDMFT (1-electron density 
matrix, orbital phase). 

2. A general application of TD-DFT to photochemistry will require more so- 
phisticated functionals, depending probably on additional parameters other 
than just the density. Studies of TD-RDMFT indicate that TD-DFT func- 
tionals may also need as variables the occupation number and the time- 
dependent phase of the orbitals. 
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Figure 1: Schematic representation of potential energy surfaces for photophysical 
and photochemical processes: Sq, ground singlet state; Si, lowest excited singlet 
state; Ti, lowest triplet state; ABS, absorption; FLUO, fluorescence; PHOS, phos- 
phorescence; ISX, intersystem crossing; CX, conical intersection; TS, transition 
state. 




Figure 2: Schematic of excitations in the two-orbital two-electron model. Exci- 
tations on the left hand side of the diagram are spin-flip excitations while those 
on the right hand side are spin-preserving excitations. 
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Table 1: Summary of formulae within the two-orbital two-electron model 
(Fig. [2]). /i is minus the ionization potential of orbital i. Aa is minus the 
electron affinity of orbital a. Aa{i~^) is minus the electron affinity of orbital a 
for the ion formed by removing an electron from orbital i. ujt, ^S: and com are, 
respectively, i ^ a excitation energies to the triplet, singlet, and mixed symme- 
try states. The mixed symmetry state has energy ojm = Aa{i~^) — h- ASCF 
quantities are obtained by the usual multiplet sum procedure (|12p except that 
a truncated Taylor expansion of the xc-functional has been used in the DFT 
case ()13p . The identification of /j and Aa in the TD-DFT case is based upon 
OFF theory 



ASCF HF 

U = ei 

Aa = Ea 

Aa{i~^) = Aa - {aa\fH\n) + {ai\fH\ia) 
= ea-ei- {aa\fH\n) + {ai\fH\ia) 

UT = LOM - {ia\ fH\ai) 

oJs = + {ia\fH\ai) 

hnearized ASCF DFT 

Ii = ei- \ {n\fH + flc\n) 

Aa = ea + \{aa\fH + fl'^\aa) 

Aa{i-^) = Aa- {aa\fH + fp^\n) 

= ea-ei + i(aa - ii\fH + /li^|oo - n) 

UT = UJM + {aa\f^;^ - f^;}\n) 

us = UJM - {aa\ftc - fPJ'lii) 



CIS (TDA-TD-HF) 

li = ei 

Aa = Ea 

Aa{i'^) = Aa - {aa\fH\n) + {ai\fH\ia) 
= ea-ei - {aa\fH\n) + {ai\fH\ia) 
ujT = UM - {ia\fH\ai) 
(^s = + {ia\fH\ai) 

TDA-TD-DFT 

It = ei 

Aa = ea + {aa\fH\n) + (ai|/2^|m) 

Aaii"^) = Aa - {aa\fH\n) + (ai|///|m) 

UJM = ea- ei + {ai\fH + /li^|«a) 

ut = UJM - {M fu + fl'^\ai) 

ujs = UJM + {ia\fH + ftc\ai) 
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Figure 3: The front (black) ladder of Jacob's jungle gym ([22|) is Jacob's ladder 
for Exc (|23p . Spin indices have been dropped for simplicity, however present 
practice is to always use spin-density functional theory where the density is a 
two-component object, {pa^pp)- LDA, local density approximation, involves only 
the density, p(r), at each point; GGA, generalized gradient approximation, in- 
volves both p{y) and the reduced gradient, x{r) = |V/9(r)|//?^/^(r); mGGA, meta 
generalized gradient approximation, involves p(r), x(r), and the local kinetic en- 
ergy, T(r) = 'T-p^p(r)V^'0p(r), or the Laplacian of the density, V^/o(r); in 
climbing to the fourth rung — that of hybrid functionals, exact exchange, and 
related functionals — explicit information about occupied orbitals is also incor- 
porated into the functional; finally on the fifth and highest rung — that of double 
hybrids, functionals based upon the adiabatic connection fluctuation-dissipation 
theorem, and related functionals — explicit information is added about unoc- 
cupied orbitals and their orbital energies. Effectively at the highest level DFT 
closely approaches conventional many-body theory. 
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Figure 4: Schematic of spin-flip excitations in the two-orbital two-electron model. 



